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^ : Abstract 

00 
O 

' We study the energy dependence of the cross sections for nucleon removal of 

projectiles. It is shown that the Glauber model calculations with nucleon- 

Q ' nucleon t-matrix reproduce well the energy dependence of the breakup cross 

^ ' sections of ^B. A DWBA model for the breakup cross section is also proposed 

: ' and results are compared with those of the Glauber model. We show that to 

rS 

■ obtain an agreement between the DWBA calculations, the Glauber formalism, 

and the experimental data, it is necessary to modify the energy behavior of the 
effective interaction. In particular, the breakup potential has a quite different 
energy dependence than the strong absorption potential. 
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I. INTRODUCTION 



Study of breakup reactions with halo nuclei is one of the main tools to understand their 
structure. The measurements of the width of momentum distributions of fragments, the 
magnitude of the total reaction cross sections, and single- and double-nucleon removal cross 
sections have been of major usefulness to unravel their internal properties (for a review, see, 
e.g., 1^). These measurements have been carried out at relatively high energies, in the beam 
energy range of 30 MeV/nucleon to 1.2 MeV/nucleon. 

The Glauber formalism is the major theoretical approach in use to analyse these mea- 
surements. This formalism is well established and yields very reasonable results for the 
reactions involving stable nuclei at high energies. In particular, a direct connection of the 
quantum mechanical breakup amplitudes and semiclassical calculations can be done in the 
Glauber formalism in a very intuitive way . 

In perturbation theory the transition amplitude is given by 

T;, = (0^vl/-|t/|0,vl;+) , (1.1) 

where < (< (f)f\) denotes the initial (final) internal wavefunction of the nuclei, \I'~ (^^) 
is the incoming (outgoing) scattering wave of the center of mass, and U is the interaction 
potential. The Glauber formalism uses eikonal wavefunctions for the scattering waves. The 
product vl/^*.vl/+ is then simply a plane wave displaced by a (eikonal) phase which is directly 
proportional to the integral of the absorptive potential along the beam direction, the z-axis. 

The use of eikonal wavefunctions is a crucial step in the Glauber formalism. Indeed, as 
shown by Glauber in his excellent lecture notes 0], the z-integration can be done by parts, 
and the potential U in eq. (1) will only appear in the exponential phase. This procedure 
is valid if the excitation potential U is the same as the absorptive potential appearing in 
the eikonal phases. This is the case for example in the calculation of total reaction cross 
sections. However, it is not a general premise. For example, in inelastic excitations of surface 
vibrations, the absorptive potential, Uabs, causing the phase-shifts in the elastic channel is 
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not the same as the interaction potential Uex which leads to surface vibrations, although 
(in some models) they can be related by derivatives. But, even in such situations, the 
energy dependence of both potentials are roughly the same. However, while the excitation 
potential is related to a few reaction channels, the absorptive potential carries information 
of all channels which may lead to the absorption of the scattering waves. Thus, one expects 
that a difference in the energy dependence of the interaction and the absorptive should be 
manifest in some sensitive cases. 

A good place to look for a deviation from the Glauber theory is the breakup reactions 
involving halo nuclei. This is because the energies involved in the breakup are basically the 
separation energies of the valence nucleons while the core nucleons which are also relevant 
for the absorptive part of the potential have much larger separation energies. Also, the 
spatial distribution of the valence and core nucleons are very different so that they influence 
differently on the absorptive (for which all nucleons participate) and the excitation (for 
which only the valence nucleons participate) potential. 

In section 2 we illustrate the connection of the Glauber formalism and the nucleon- 
nucleon cross sections. We apply the theory to the calculation of proton removal cross 
sections from projectiles at several bombarding energies. In section 3 we develop a 
DWBA formalism for the stripping reactions. In section 4 we show how to relate the DWBA 
calculation with the Glauber formalism, using effective interactions. Our conclusions are 
presented in section 5. 

II. GLAUBER MODEL AND NUCLEON-NUCLEON SCATTERING 

AMPLITUDES 

In the Glauber theory, after the z- integration, the remaining integrals in eq. (1) can be 
easily related to the concept of impact parameter and to absorption and survival probabil- 
ities. For example, simple manipulations show that the nucleon removal cross sections in 
high energy collisions are described in the Glauber theory by 
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a — 27: J dbb [1 — exp {—2TmXv)] exp (— 2TmXc) , (2.1) 

where Im stands for the imaginary part, v (c) denotes for valence (core) particles, and x 
are the eikonal phases given by 

X^ib) = -^ J^Jz'Uir'), (2.2) 

where v is the projectile velocity, r' = vP+^, and U is the optical potential for the system 
composed of the particle i {— v, c) and the target. The term inside the brackets in eq. (2) 
can be interpreted as the probability that the valence nucleons will be removed in a collision 
with impact parameter 6, while the exponential term outside brackets is the probability that 
the core nucleons will survive. This product integrated over all impact parameters gives the 
cross section for (valence) nucleon removal. 

A great simplification introduced by Glauber was to relate the optical potentials to the 
nucleon-nucleon cross sections. This can be done easily by noticing that the optical theorem 
for the forward nucleon-nucleon amplitude yields 

W(q = 0) = f{e = 0°) = a^,,{E) [1 - ia{E)] , (2.3) 

where (Tjv at is the nucleon-nucleon cross section and a is the real-to-imaginary ratio of the 
nucleon-nucleon scattering amplitude. With the assumptions that only very forward angles 
are involved, and that the nucleon-nucleon interaction is of very short range (i.e., a delta 
function interaction) , one can construct optical potentials for the nuclear scattering in terms 
of the folding integrals 

Ui{R) = t7Viv(q = 0) I p,(r)p^(R - r) £r , (2.4) 

where pi and pa are the ground state densities of the projectile i, and the target A, respec- 
tively. 

As an application of this model, let us consider the proton removal cross sections of 
projectiles in reactions with carbon targets. For the valence nucleon we get the density 
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distribution from a Woods-Saxon+spin-orbit potential model for a proton in the J93/2 orbital 
(the parameters are given in section 4). For the core CBe) density we use the ground state 
density parameterized as p(r) = po[l + cr"^ /a?] exp (r^/a^), with a = 1.77 fm and c = 0.327 
fm. The result of the calculation is shown by the short-dashed line of figure 1. The ex- 
perimental data are from ref. [Q. Although the magnitude of the cross section is a little 
overestimated, we see that the energy dependence follows very closely that of the experi- 
mental data. The dashed curve is the calculation renormalized to the lowest experimental 
data point. In fact, the reasonably good agreement between the energy dependence deduced 
from the Glauber theory and the experimental data on the total nuclear cross sections, and 
nucleon removal cross sections, is well established, both for stable and unstable nuclei. 

Also shown in figure 1 (solid curve) is the model developed by Hansen In his model, 
the nucleon removal cross section is forced to have the same energy dependence as the 
total reaction cross section. The total reaction cross section has a slightly different energy 
dependence than the valence nucleon removal cross section. This can be best seen from the 
Glauber theory. The calculation of the total reaction cross section amounts in replacing 
the integrand in eq. (2) by 6 [1 — exp {—2XmXaA)] where now XaA is the eikonal phase for 
the collision of the projectile a and the target A. In fact, we see that the Hansen's model 
predicts a rather different energy dependence of the proton removal cross section. The data 
favor the calculation following eq. (2). 

From a general point of view, the energy dependence of the total and the nucleon removal 
cross sections are directly related to the underlying optical potentials for the reaction. From 
the above discussion we can see that these optical potentials should have a similar energy 
dependence as the nucleon-nucleon cross section. To study this idea further, let us formulate 
a DWBA model for the breakup cross section. The use of an effective nucleon-nucleon 
interaction, M3Y interaction, will serve as guide to understand the link between the optical 
potentials and the nucleon-nucleon cross sections. 
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III. DWBA BREAKUP AMPLITUDES 



Let us consider the general case of the stripping of the projectile a incident on target A: 

a + A — >b{a-x) + B (A + x) . (3.1) 
The hamiltonian for the system is 

H = Ha + HA + TaA + VaA = H, + Hb + T.b + , (3.2) 
The transition matrix element for this reaction is given by 

T = (^"(kb, rbs)(/)b(Sb)0B(C, ^x, ^xa) \Ubu\ *a (ka, raA)0a(S6, s^, r6^)0A(O) > (3-3) 

where 0a, 0b, 0a, <t>B are eigenstates of Ha, H^, Ha, Hb, respectively, i.e., Ha(Pa — ^A(f>A, 
a — ^a(pai ctc, ^ are distorted waves of the particles a and b, i.e., {TaA + UaA)^a — 
{E — ea — e^)^^. The internal coordinates of a, b and x, respectively, are denoted by Sj 
{i — a, b, x), Vij are the relative coordinates of particles i and j, and ^ is the internal 
coordinate of particle A. 

We use the coordinate relationship 

Tbs = + ( — ) v^A , r^A = r^A + {—] , (3.4) 
and we integrate over the internal coordinate, ^, oi A 

i^xisx^r^A) ^ j (t>B{i, s^,r^A) 0a(O , (3-5) 
and over the internal coordinates of b and x, 

j dSb dSx 06(S6) Ipx (Sx, r^A) 0a(S6, S^;, r^^) = Cbx 0a(r6x) ^xi^xA) , (3.6) 

where \Cbx\'^ is the spectroscopic factor. 
We get for the transition matrix element 

T^Cbx J d^nx d^VbA "^b* (kb, + ~^^^) ^x{rxA) 

X Ubu{r^A, rbx) (j)airbx)'^a (K,rxA + —rbx] ■ (3.7) 



The potential Uex for the breakup channel is given by 



= Ul,A{rxA + nx) + UxAi^xA) - UaA ^xA H 

Vila 



mi, 



(3.8) 



If we now integrate over r(,^, we can use the fact that the bound state wavefunction is 
peaked at small r^a. values, so that 

/ d^nx '^l* ( kb, r;,^ + — r^^ ) U^^ijc^A, rbx)(pa{rbx)'^a ( ka, rxA + —Ybx] 
J V niB J \ rua J 

fa (kb, -^^xA^ (ka, ^xa) J d^Hx Ubu{rxA, nx)Mnx) (3.9) 
We now define a "transition" , or "excitation" , potential as 

Uexi^xA) = j d^Vbx M^bx) Ubu{rxA: ^x) , (3.10) 

SO that 

T^Cbx I dhxA (kfe, ^r,^) *-(r.A) U,x{Vxa) ^+ (k„, r,^) . (3.11) 

The above equation is our main result. It gives the t-matrix in terms of the scattering 
waves of particle a, 6, and x, and a "transition" potential Uex- This potential contains the 
information on the structure of particle a. 

If we are only interested in particle 6, assuming that the particle x is not observed, we 
can use the closure relation 



E *^^*(r.A) ^f^\v'xA) = 5 {r^A - v'xa) , 



(3.12) 



to obtain 



E|T|' = |a.|' fd^TxA 



( k„ ^TxA 

rriB 



Uex{rxA)\^ ^'a {K, ^xa) 



(3.13) 



The calculation becomes very transparent if we use eikonal functions for the distorted waves: 



cxp 



.rriA. , . ^a, 

TUB yrriB 



*a = exp [ik^.r^.^ + ixa (bxA)] , 



(3.14) 



with the eikonal phases given by 



Xb 0.x 

Xa {Ka 



1 

1 fZxA 



tiVb Jz^A yrriB 



xA ) ""^xA 



hVa J~ 



UaA {r'xA) dz'^A 1 



(3.15) 



where r^^ — \f^A^rz^^. Note that, since only the moduh of the \I/'s enter into eq. (18), 
we get 



^ \T\^ = ICfe^l^ j d^r^A Sa{rxA) Sbir^A) \Uex{rxA) 



where 



Sb{rxA) = exp 
Sa{rxA) = exp 



J —oo 



hVa J-oo 



xA 



FinaUy, the cross section for removal of particle x from a is given by 



a — 



rriamb kb Ek. 



X, spins 



(3.16) 



(3.17) 



(3.18) 



TTh^ K (2J^ + 1)(2J„ + 1) ' 

and we shall assume that Vb ~ fa, valid for high energy collisions and small binding energies 
of the incident projectile. 



IV. WAVEFUNCTIONS AND EFFECTIVE INTERACTIONS 

The ground state wavefunction of ^5, in a given magnetic substate, M, is taken as 

0f Hrft.) = E {3mlAMA\JM) ct>,^{v) \IaMa) , (4.1) 

m,MA 

where \IaMa) is the wavefunction of the ^ Be {I a = 3/2~), and 0^-,^ is the single-particle 
wavefunction of the proton j"^ — 3/2~, coupled to a total angular momentum = 2"*". 
Thus, the potential Uex in eqs. (15-17,21) depends on the initial orientation of and the 
target, what means that eq. (23) carries an average over the magnetic substates of these 
nuclei. 
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Using the properties of the Clebsch-Gordan coefficients, and the orthogonahty of the 
core wavefunctions, we get for the spin averaged potential 

= (01 E \Ue.{rxA)? = ^ / d'n.^^ UUr.A, n,) , (4.2) 

I^ZJa + -LJl^-'a + -Lj spins •' ^bx 

where Rj{r)/r is the radial part of the single-particle wavefunction (pjrn- The cross section 



IS 



2ml 



l^xf / dhhSab{h)J',.,{h) , (4.3) 



where Sah = SaSb (we neglect the small dependence of Sab on z) and 

^ex(&) = r dz iVW+^) . (4.4) 

Now we need to determine the optical potentials to proceed with the calculation. Usually 
these optical potentials are obtained from elastic scattering experiments. But, for unstable 
nuclei the situation is quite different. One generally has to construct these optical potentials 
theoretically from effective nucleon-nucleon interactions. Among these, one of the most 
popular is the M3Y interaction, which has been shown to work quite reasonably for elastic 
and inelastic scattering of heavy ions at low and intermediate energy collisions . 

In its simplest form the M3Y interaction is given by two direct terms with different 
ranges, and an exchange term represented by a delta interaction: 

t{s)=A-— + B-— + C5{s), (4.5) 

where A=7999 MeV, B = -2134 MeV, C = -276 MeV fm\ /?i = 4 fm-\ and p2 = 2.5 
fm~^. The real part of the optical potential is obtained from a folding of this interaction 
with the ground state densities of the nuclei: 

Ui,{K) = J rfVi d\2 PA{ri)pj{r2) t{s) , (4.6) 

with s = R + r2 — ri. The imaginary part of the optical potential is usually parameterized 
to be JmU = XUmsy, with A = 0.6 — 0.8. 
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The M3Y interaction (28) has been modified to account for the energy dependence on 
the beam energy. However, for the energy range of 5 — 50 MeV/nucleon, only a small energy 
dependence was introduced |0 as a variation of the exchange term. 

To study the breakup of projectiles, we will use the form given by eq. (28) for 
the M3Y interaction with the ^Be and ^B densities as in section 1, and a proton gaussian 
density of radius equal to 0.7 fm. The radial wavefunction of the proton, -R3/2, was obtained 
in a Woods-Saxon+spin-orbit potential, i.e., V{r) = Vo[l — Fso(}-s){ro/r)d/dr] f{r), with 
/(r) = [1 + exp ((r — R)/a)]^'^ with parameters Vq = —44.66 MeV, a = 0.56 fm, ro = 1.25 
fm, R = 2.391 fm, Fgo = 0.351 fm, which reproduces the binding energy, e = 0.138 MeV, of 
^B. The spectroscopic factor is taken as unity. 

In figure 2 we plot the function jFg^(6), which contains the information not only of the 
ground-state wavefunction of the ^B, but also on the effective interaction. In figure 3 we 
plot the profile function Sab{b), which depends only on the effective interaction. We have 
calculated it for the energies E/A =30, 150, 300, 800, and 1200 MeV, respectively. The 
magnitude of the cross section is proportional the area below Sab{b) x J-'ex{b). Since J^ex does 
not depend on the beam energy, the energy dependence is solely due to Sab- Since the M3Y 
interaction does not depend on the energy, the energy dependence is a consequence of the hv 
factors in the denominators appearing in eqs. (22). This causes the nuclear transparency, 
described by the factor Sab, to increase for small b as the beam energy increases. As a 
consequence, the cross sections increase with energy. A comparison with the experimental 
data in figure 4 (dotted curve) shows the departure of the calculated cross sections from the 
experimental data at large energies. 

It is clear that we have to modify the effective interaction in eq. (29) so as to incorporate 
the energy dependence. A simple way to do that is to make (29) have the same energy 
dependence as in eq. (4,5). It should be noticed that the potential Uex in eq. (16) is not 
the same as the potentials appearing in the phase of the scattering waves, in the sense 
that it does have neither the same magnitude, nor the same spatial dependence. If we take, 
although it is not necessary, the same energy dependence as in (4) we can obtain the effective 
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interaction as 

fiv 

t(E, s) = -i— aMN{E) [1 - ia(E)] t(s) , (4.7) 

where to — 421 MeV is the volume integral of the M3Y interaction. Note that eq. (30) 
gives the same removal cross section as the M3Y interaction for E==30 MeV. Inserting this 
result (eq. 30) in eq. (29), we can determine the imaginary part of the optical potential 
automatically. 

We repeat the calculation for the proton removal cross sections of using the effective 
interaction (30) in eq. (29) and the calculated cross sections by eq. (26) are shown in figure 
4 (dashed hue) . We see that the energy dependence of the cross section changes drastically 
and follows more closely the trend of the experimental data. 

As mentioned above we do not need to assume that the absorption potential and the 
excitation potential have the same energy dependence. It is reasonable to assume that the 
absorption potential follows the receipt of eq. (29), with t given by (30), since this has the 
same energy dependence occurring in all calculations based on the Glauber formalism for 
total reaction cross sections, which are known to agree reasonably with the experimental 
data. Thus, we change the excitation potential to adjust its energy dependence to the data 
points. We find out that a simple energy dependence of the form t{E,s) oc E'^'"^^ t{s), 
for E < 200 MeV, and Z4x constant .t{s), for E > 200 MeV, reproduces the trend of 
the experimental data, as we show in figure 4 by the solid line, with a normalization factor 
which best fits the data. 



A. Conclusions 

In summary, it is found that the energy dependence of the experimental removal cross 
sections of can be obtained by Glauber model calculations with the nucleon-nucleon 
t-matrix. We have studied also the relation between the Glauber model and the DWBA 
formalism. The DWBA approach to nucleon removal cross sections in general agrees with 
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Glauber calculations if the optical potential in the elastic channel has the same energy 
dependence as the breakup potential. For halo nuclei this is not necessarily true, as we have 
shown for the breakup of the nucleus. This finding might have important consequences, 
not only for the breakup of halo nuclei but also for their excitations to bound states. More 
studies with halo nuclei are needed in order to clarify the role of effective interactions in 
the construction of optical potentials, and of their connection to nucleon-nucleon scattering 
amplitudes. 
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Figure Caption 

Fig. 1 - Proton removal cross sections of projectiles on carbon targets as a function 
of the incident beam energy. Data points are from ref. 0. Solid line is a calculation based 
on a model by Hansen 0. Short-dashed line is a calculation based on the Glauber model, 
eq. (2). Dashed line is simply the short-dashed line downshifted by a factor of 0.83. 

Fig. 2 - Excitation function J-'^x (see definition in text; eq. (27)), as a function of the 
impact parameter. 

Fig. 3 - Nuclear transparency function for the proton removal of projectiles at several 
energies incident on carbon targets. 

Fig. 4 - Proton removal cross sections of ^B projectiles incident on carbon targets as a 
function of the incident energy. Dotted line is the result of a DWBA calculation with the 
folding potentials with M3Y interaction. For the dashed line the effective interaction was 
taken with the same energy dependence as the nucleon-nucleon scattering amplitude. Solid 
line is the result obtained with a fit for the energy dependence of the breakup potential, 
different from that of the absorption potential. For details, see the text. 
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